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We present a first-principles theory of dynamical spin excitations in the presence of spin-orbit 
coupling. The broken global spin rotational invariance leads to a new sum rule. We explore the 
competition between the magnetic anisotropy energy and the external magnetic field, as well as the 
role of electron-hole excitations, through calculations for 3d-metal adatoms on the Cu(lll) surface. 
The spin excitation resonance energy and lifetime display non-trivial behavior, establishing the 
strong impact of relativistic effects. We legitimate the use of the Landau-Lifshitz-Gilbert equation 
down to the atomic limit, but with parameters that differ from a stationary theory. 
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I. INTRODUCTION 

The understanding and design of new technologies 
based on magnetic materials, in the fields of spintron- 
ics and magnonics, begets and profits from quantita¬ 
tive theoretical approaches. Recently, the central role 
played by spin-orbit coupling (SOC) has been serendipi- 
tously revealed. It underlies many concepts for the gen¬ 
eration, manipulation and detection of spin currents^, 
also at the nanoscale: spin dynamics, current-induced 
magnetization switching and magnetic stability^^"*. The 
energy gap in the spin excitation spectrum, from ex¬ 
tended materials down to a single adatom, is the result 
of SOC, being measurable, for instance, with ferromag¬ 
netic resonance (FMR)® or inelastic scanning tunneling 
spectroscopy (ISTS)®^^^. Physical intuition suggests a 
connection between the gap and the magnetic anisotropy 
energy (MAE), while application of a dc external mag¬ 
netic field (R™*) should lead to a gap change of 
(Zeeman shift). Stoner (electron-hole) excitations, to¬ 
gether with SOC, contribute to the spin dynamics, their 
impact on the energy gap, Zeeman shift and excitation 
lifetime being expected but fairly unexplored. 

Up to now, all ah initio studies on transverse dynam¬ 
ical spin excitations have neglected SOC. Furthermore, 
the addressed magnetic states are collinear {e.g. ferro¬ 
magnetic), since the complexity, both theoretical and 
computational, increases dramatically in the general 
case: charge excitations, longitudinal and transverse spin 
excitations may couple in a non-trivial manner (see dis¬ 
cussion in Ref. 12). From the tight-binding perspec¬ 
tive, the work of Costa et al.^^ has proven invaluable 
in characterizing the impact of SOC on the spin-wave 
dispersion and lifetime, going beyond the adiabatic ap¬ 
proximation^^. Two kinds of theoretical approaches 
build upon first-principles electronic structure calcula¬ 
tions. One type is based on many-body perturbation 
theory (MBPT)^^^^^, constructing the non-interacting 
Green function (GF) from the density functional theory 
(DFT) eigenstates. The other type is based on time- 
dependent DFT (TDDFT)^®d9^ lias been applied to 


bulk systems^®’^^’^*^“^^, thin films^"^ and adatoms on sur- 
faces^^“^^, with pioneering work on SOC^®, but still not 
starting from a spin-polarized ground state^®. Recently 
we provided a connection between MBPT and TDDFT, 
to describe the interaction between electrons and spin 
excitations®*^; this link was also found for the theory of 
spin-fluctuation-mediated superconductivity®®. 

Here we present a method for the calculation of dy¬ 
namical magnetic response functions based on TDDFT 
and incorporating SOC. This new scheme is imple¬ 
mented within the Korringa-Kohn-Rostoker (KKR) GF 
approach®®. Of interest is the ability to treat an external 
magnetic held (R®’^®) in any direction, to investigate non¬ 
trivial orientations of the magnetic moments and the role 
of the MAE. We derive a novel magnetization sum rule 
that constrains the exchange and correlation (xc) kernel 
when SOC is present, which is essential for the theory and 
the calculations. By investigating 3d metal adatoms on 
the Cu(lll) surface, with a focus on the experimentally 
studied Fe adatom®, we demonstrate that the gap in the 
excitation spectrum is connected to the MAE through 
spin dynamics parameters, with electron-hole excitations 
taking center stage. We also show that, varying the ori¬ 
entation of R®’'®, the properties of the spin-excitations, 
i.e. excitation energies, lifetimes and (/-factors, are non¬ 
trivial and anisotropic, depending on the orientation of 
the magnetic moment. 

Our paper is organized as follows. In Sec. II the for¬ 
malism for the dynamical magnetic susceptibility is pre¬ 
sented, as is the new magnetization sum rule that applies 
for the SOC case. Sec. Ill describes the computational 
details. The results for 3d adatoms on Cu(lll) are dis¬ 
cussed in Sec. IV, and the interplay between MAE and 
SOC is explained in Sec. V. Our conclusions are given in 
Sec. VI, while three appendices provide further details. 


II. DYNAMICAL SUSCEPTIBILITY 

We begin by seeking the change of the spin density 
matrix, dp = dnao + Srh ■ a, due to an external time- 
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dependent potential of the same form, = 6Vq^^ (To+ 

SB°^^ ■ a. In linear response, 

Sp{r;t) = Jdf'Jdt'x{r,f';t-t') ■ 6V^^^ir']t') , ( 1 ) 

with a the vector of Pauli matrices and cto the 2x2 unit 
matrix. Within TDDFT, after Fourier transforming in 
time, the full response function, Xy is given in terms of 
the Kohn-Sham (KS) response, fh® Hartree- 

exchange-correlation (Hxc) kernel, through 

X(a;) = [1 - x^H^) x^H^) (2) 

(here and in the following spatial dependence and inte¬ 
grations are omitted for brevity). We add SOC^^ and 
spin and orbital Zeeman couplings to an external static 
field to the KS Hamiltonian, self-consistently, 

•HKS .^ext _ (3) 

The KS response is written in terms of the KS GF, 
G(E) = {E-n^^y\ as 

1 pEf 

Xa|(w) = TT[aa G{E + OJ + iO) crplmGiE) 

+ a a Im G{E) ap G{E — uj — iO)] , (4) 

with a,fi = x,y,z. The trace is over the spin compo¬ 
nents, and 2iImG(F;) = G(F;-kiO)-G(£;-iO). 
is the sum of the Hartree kernel, 2/|r — r'|, and of the 
xc kernel, which contains the many-body effects. In 
the adiabatic local spin density approximation (ALSDA), 
adopted in this work, the xc kernel is local in space and 
frequency-independent, being a function(al) of the parti¬ 
cle and spin densities only. 

In the absence of SOC, the full response function de¬ 
couples into a transverse and a longitudinal part. The 
former describes damped precessional motions of the 
magnetization, while excitations that change the mag¬ 
nitude of the charge or spin densities are in the lat¬ 
ter. Global SU(2) invariance implies the existence of a 
Goldstone mode: the spin density is an eigenfunction of 
X“^(tLi = 0) with vanishing eigenvalue. This is of utmost 
importance in numerical calculations, as small inaccura¬ 
cies in the KS susceptibility and in the xc kernel shift 
the Goldstone mode to a finite frequency, in the meV 
range^^’^®, where the gap opened by SOG is expected. 
Gorrective schemes were proposed^^’^®’^® by adjusting 
one or both of these quantities to place the Goldstone 
mode at zero frequency, when SOC is neglected. Next 
we derive a new sum rule, connecting the spin density to 
the xc magnetic field, SOC and external magnetic field, 
so that the gap arises unambigously from the latter two. 

Suppose that the xc magnetic field, lies in the 
0 -direction, defining the local spin frame of reference®'*’. 
Then the spin density is given by (see Appendix A): 

1 

m, = --lmTT dEG^yE)A{E)Gu{E) . ( 5 ) 


The effective spin splitting is, cf. Eq. (3) and Appendix A: 

where the auxiliary GFs are Ga{E) = (A — 
with cr = tj4-- Identifying the spin-flip KS susceptibil¬ 
ity via Eq. (4), X+_ = [xxx - ^Xxy + '^Xyx + Xyy)/‘^y 
we can relate the spin density to the xc magnetic field, 
= 2 x+® (w = 0) B^'^ + Srriz , where Sniz arises from 
all contributions to A, excluding the xc part. The trans¬ 
verse xc kernel in the ALSDA is just K± = 2B^'^/mz, 
yielding the magnetization sum rule 

mz{f) - Jdr' x+-ir,r';0) Kj_{f')mz{f') = (5m^(f) . 

(7) 

When Smz{r) = 0 (no SOC or external field), the de¬ 
nominator of Eq. (2) vanishes and the Goldstone mode 
is recovered. If only an external field is applied, the spin 
excitation is located at w ~ 2Bf^^ (Zeeman shift), while 
SOC gives rise to a gap even in zero field. If Eq. (7) is 
not satisfied due to numerical inaccuracies, K±(f) is ad¬ 
justed such that the sum rule holds, using as input the 
calculated 6mz{r), rrizif) and x+-(7’j = 0). 

Next we briefly summarize the computational details, 
before presenting results of our TDDFT formalism for 
adatoms on Cu(lll). 

III. COMPUTATIONAL DETAILS 

The GFs are evaluated by the KKR method®^ where 
the electronic structure of the adatoms is computed in 
two steps. First the electronic structure of a 22 layer Cu 
slab is calculated. Then each adatom is self-consistently 
embedded on the surface of this slab, in real space, to¬ 
gether with 12 nearest-neighbor sites. We consider the 
Fe adatom relaxed to the surface by 10% of the interlayer 
distance of Cu and use, for the sake of comparison, the 
same vertical distance for all adatoms. We employ the 
rigid spin approximation, whereby the direction of the 
spin density is taken to be collinear inside each atomic 
sphere®®. To tackle canting of the spin moment due 
to perpendicular external magnetic field and anisotropy 
easy axis, its direction is updated during the iterations 
until self-consistency is achieved. The MAE is calculated 
by band energy differences following the magnetic force 
theorem®®, 

. ( 8 ) 

For the susceptibility calculations, the proposal of Lounis 
et aZ.^®’^® is extended to an spdf basis built out of regu¬ 
lar scattering solutions evaluated at two or more energies, 
by orthogonalizing their overlap matrix (see Appendix B 
for details). This basis can reproduce ground state data 
reliably and is used to include SOG in the GFs, self- 
consistently The spatial dependence of the susceptibil¬ 
ity is restricted to the magnetic adatom; tests including 
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neighboring copper atoms show no significant impact on 
the transverse spin excitations. The energy integration 
in Eq. (4) is performed as detailed in Refs. 25 and 26. 


IV. SPIN EXCITATIONS OF Cr, Mn, Fe AND 
Co ADATOMS ON THE Cu(lll) SURFACE 


As an application of our method, we explore the spin 
excitation spectra of 3d adatoms on the Cu(lll) surface. 
The groundstate properties and the dynamical spin exci¬ 
tation spectra are described separately. 


A. Ground state properties 



E-Ep(eV) 


FIG. 1. Atom-projected total density of states for Cr, Mn, Fe 
and Co adatoms on the Cu(lll) surface (positive for majority 
and negative for minority spin). Energy measured from the 
Fermi energy of the substrate. 


The adatom-projected local density of states for Cr, 
Mn, Fe and Co adatoms on the Cu(lll) surface is shown 
in Fig. 1, while ground state properties are listed in Ta¬ 
ble I. The spin moment, m®, is maximum for Cr and Mn, 
and then decreases steadily for Fe and Co. The orbital 
moment, m°, is small for Cr and Mn, a consequence of 
the almost half-filled d-states, and large for Fe and Co, 
due to partial occupation of the minority d-states. As 
expected, Cr and Mn adatoms with nearly half-filled d- 
shells have larger spin magnetic moments (nf) and lower 
orbital magnetic moments (m°) than those of Fe and Co 
adatoms. The small magnitude of the orbital moments 
for Cr and Mn correlates with small MAEs, the converse 
being true for Fe and Co. The preferred orientation of 
the magnetic moments of Fe and Co is normal to the 
surface while Cr and Mn lie in-plane. 


B. Dynamical spin excitations 

Now we proceed to the calculations of the dynamical 
magnetic susceptibility. Eq. (2) is solved using the full 
basis expansion of the GFs, but for discussion we define 
an adatom-averaged quantity, corresponding to the net 
response to a site-dependent TD external magnetic field: 

dfj dr'Xai3{r,r';u}) . (9) 

First we discuss the impact of SOC on the adatom- 
averaged KS susceptibility. (w) is linear for small fre- 




Cr 

Mn 

Fe 

Co 

m® (pb) 

4.07 

4.31 

3.23 

1.97 

m° (^b) 

-0.02 

0.02 

0.55 

0.52 

Ea (meV) 

-0.29 

-0.33 

4.96 

2.25 


TABLE I. Spin (rrd) and orbital (m°) moments for Cr, Mn, 
Fe and Co adatoms on the Cu(lll) surface, and MAE, see 
Eq. (8). All values for the easy axis conhguration (in-plane 
for Cr and Mn, out-of-plane for Fe and Co). 



FIG. 2. Density of transverse spin excitations for Cr, Mn, Fe 
and Co adatoms on the Cu(lll) surface (R®’'* = 0). (a) from 
Imx+?.(ta), Eq. (4). (b) from ImX+_(a)), Eq. (2). Dashed 
lines show the result of the estimate based on the adiabatic 
approximation (see main text). 


quencies. Fig. 2(a) shows Imy^® (w), describing spin-flip 
excitations between occupied and empty states (Stoner 
excitations). SOC is found to have a negligible impact 
on (w) for the Cr and Mn adatoms, in line with their 
low orbital magnetic moments and MAEs. For Fe, SOC 
causes a noticeable change only on x+-(<j-’ = 0) (^0.2%), 
while for Co there is also an increase in the slope of 
Imx+® (w) by 12%. 

Consider now the full response function. The longitu¬ 
dinal and transverse parts of the magnetic susceptibility 
are weakly coupled by SOC. The imaginary part of the 
dominant eigenvalue of x(<^) corresponds to the density 
of states of the main magnetic excitation, and is shown 
in Fig. 2(b) for all adatoms. Although the resonance en¬ 
ergies (R®’'* = 0) follow the trend of the adatoms’ MAE 
(Table I), their positions are strongly shifted from what 
is expected from a simple approximation (dashed lines in 
Fig. 2(b); discussed in the following). 

We characterize the main resonance in x(<^) through a 
simple model of the spin dynamics, the Landau-Lifshitz- 
Gilbert (LLG) equation^"^, widely used for larger mag¬ 
netic systems: 


dm® 

dt 


—ym® X i?®®-|-ry 


m® dm® 
m® dt 




dE 

drn^ 


( 10 ) 


with ?fi® the spin moment, m® its length, 7 the gyromag- 
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netic ratio (equal to 2^b for a free electron) and r] the 
damping parameter. We set S = 1 and absorb /j,b in i?®®, 
so 7 = 2 . As discussed by Kittel^® and found from our 
calculations, the orbital moment mostly follows the spin 
moment, m° oc m®, so it is not considered independently. 

Suppose Cz in equilibrium, so 

and we linearize Eq. (10) adding a small time-dependent 
perturbation in the xy-plane. From Imx_|_(w), the res¬ 

onance location, Wmax) is given by (Appendix C), 


7 B 


eff 




B®® = 




2Ea 


mr 


+ 1 + ^ B®’^‘ (11) 


with a model energy of the form 
Ea 


E = -- 


r(m®-e,)"-(m® + m°)-B® 


( to ®)2 

The spectroscopic g-factor is defined as 
dcJmax 7 


9 = 




, m 

^ -^ 

TO® 


( 12 ) 


(13) 



FIG. 3. (a) Frequency cJmax at which Imy is maximum vs 

applied external magnetic held, for the Fe adatom on 

Cu(lll). The static held was applied normal to the surface, 
with spin-only (s, squares) or spin and orbital coupling (s-fo, 
diamonds); and in the surface plane, with spin-only (s, down 
triangles) or spin and orbital coupling (s-l-o, up triangles), 
(b) FWHM vs resonance energy, as extracted from the Imy 
curves (s-fo). The spin moment is normal to the surface (di¬ 
amonds), canted (down triangles) or in the surface plane (up 
triangles). Linear hts also shown. 


A non-zero orbital moment (SOC) implies an effective 
g (X j (l + ^'), which is the usual explanation for g ^ 2 
found experimentally (for instance in FMR®). Damping 
renormalizes g further by a factor 1 / -|- , so that 
a shift can be expected even without SOC (to® = 0). 
Lastly, the inverse of the full-width at half maximum 
(FWHM, r) of the spin excitation resonance can be used 
to estimate the corresponding lifetime, t « h/{2T). All 
model parameters are extracted by fitting the LLG form 
of x+_(a;) given in Eq. (C4) to the first-principles data 
in Fig. 2(b), and are collected in Table II. 

If we neglect dynamical corrections (7 = 2, ry = 0), 
the MAE leads to a resonance at 4 Bq/to® (dashed lines 
in Fig. 2(b); Ea is taken from Table I). The downward 
shift of the resonance energy is due to 7 7 ^ 2 and an ef¬ 
fective Ea smaller than the force theorem estimate (see 
Table II). 7 7 ^ 2 is determined by the electronic structure, 
as discussed in Refs. 9, 25, 26, and 39. The discrepancy 


in Ea may originate from the following: the value ob¬ 
tained from Eq. ( 8 ) is an energy difference between two 
orthogonal magnetic orientations, while the LLG model, 
see Eq. (10), suggests the dynamics depend on the varia¬ 
tion of the energy around the equilibrium direction of the 
adatom spin. An alternative explanation would be the 
large DOS peak for Fe and Co at the Fermi energy: the 
MAE computed from Eq. ( 8 ) is sensitive to peak shifts 
between the two orientations, which may arise from the 
frozen potential approximation. The damping parameter 
ry plays a minor role for Cr and Mn adatoms, but reduces 
g for Fe and Co by 5% and 22%, respectively. The spin 
excitation lifetimes are much shorter for the Fe and Co 
adatoms than for the Cr and Mn adatoms. This is due 
to the r oc Wmax Scaling of the FWHM, and to the con¬ 
nection between rj and the slope of Imx^® (w), Fig. 2(a). 
SOC enhances rj only for the Co adatom, for the same 
reason (see discussion of Fig. 2(a) in the main text). 



Cr 

Mn 

Fe 

Co 

Ea (meV) 

0.14 

0.02 

4.51 

1.37 

7 

1.63 

1.74 

1.73 

2.36 

V 

0.07 

0.05 

0.33 

0.80 

g (spin only) 

1.63 

1.74 

1.64 

1.84 

g (spin-l-orb) 

1.62 

1.75 

1.92 

2.33 

T (ps) 

22 

1200 

0.11 

0.091 


TABLE 11. Spin dynamics parameters for Cr, Mn, Fe and 
Co adatoms on the Cu(lll) surface, obtained by fitting the 
transverse susceptibility, Eq. (9), to the LLG model. Eq. (11) 
defines Ea- The yy-factor, Eq. (13), distinguishes or 

(m® + m°) ■ B®’'* couplings, r for Mn is hard to extract, due 
to a very narrow spin excitation peak. 


V. INTERPLAY BETWEEN MAGNETIC 
ANISOTROPY AND MAGNETIC FIELD 

Next we focus on the effect of B®’^* on the spin excita¬ 
tions of the Fe adatom. First we apply the field parallel to 
the MAE easy axis, resulting in the linear Zeeman shift. 
Wmax is shown in Fig. 3(a) for the spin-only (squares) and 
spin-borbital (diamonds) Zeeman couplings, see Eq. (3). 
In the spin-only case, g = 1.70 (1.64 using the LLG 
model. Table H). This does not depend on SOC and is re¬ 
lated to details of the electronic structure®. When the or¬ 
bital coupling is included, g rises from 1.70 to 1.97, close 
to the factor (l -b ^), from Eq. (11). Using ^ = 0.17 
and the data in Table II, yy = 1.92 from the LLG model, in 
good agreement with the first-principles data. In Ref. 8 , 
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this system was studied by ISTS. The measured gap (~ 1 
meV) is lower than the one found in our calculations 
(^4.6 meV), which may indicate we overestimated the 
MAE. The experimental g value (^2.1) is quite close to 
the one computed including spin and orbital Zeeman cou¬ 
plings {g = 1.97). The linewidth is proportional to the 
resonance energy, Fig. 3(b) (diamonds), but the slope is 
larger in the experimental data. 

When the field is applied perpendicular to the MAE 
easy axis, the equilibrium direction of the Fe spin moment 
progressively cants away from the surface normal, becom¬ 
ing parallel to (when Wmax ~ 0). Fig. 3(a) shows 
that w„iax slowly decreases for small dipping near 

the critical field, beyond which it increases again, recov¬ 
ering the linear dependence on i?®’®* for Wmax > 5 meV. 
In this regime, g = 1.81 (spin Zeeman) ot g = 1.93 (full 
Zeeman) differ from the values obtained in the previous 
linear case. This arises from the different values of m°, of 
7 and of rj, making the spin dynamics anisotropic. The 
FWHM is still linear in Wmax, Fig. 3(b) (triangles), but it 
does not extrapolate to zero as in the out-of-plane case. 
For the same Wmax, the lifetime of the spin excitation 
strongly depends on the orientation of the spin moment. 
From Fig. 3(b), with w^ax 4.6 meV, t « 110 fs (out-of- 
plane) and r « 85 fs (in-plane), which is a 20% change. 


VI. CONCLUSIONS 

In this paper we presented a detailed first-principles 
analysis of the spin dynamics of magnetic adatoms, made 
possible by an extension of TDDFT to spin-polarized sys¬ 
tems including SOC. We found a novel and invaluable 
sum rule connecting the spin density to the xc splitting 
in the presence of SOC. The key spin dynamics param¬ 
eters have been extracted from the dynamical magnetic 
susceptibility, including SOC and spin and orbital Zee- 
man terms, after mapping to the LLG model, thus legiti¬ 
mating its use down to the atomic limit. Deviations from 
standard assumptions in spin dynamics models have been 
found for the MAE (different from the one computed by 
the force theorem), the gyromagnetic ratio (7 2 , which 

may indicate spin pumping"^*^), and the origin and role 
of the Gilbert damping, rj (dominated by Stoner excita¬ 
tions, not by SOC). We also find a non-trivial behavior 
for g and the spin excitation lifetime upon application 
of a magnetic field along the easy and hard axes of the 
magnetic anisotropy. The anisotropic nature of the spin 
dynamics was established, arising from SOC and i?®^*. 
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Appendix A: Derivation of the magnetization sum 
rule for a non-spin-diagonal Hamiltonian 

We start from the Kohn-Sham Hamiltonian as given 
in Eq. (3) of the main text, 

HK®(f) + , (Al) 

and the KS Green function, G{E) = (£1 — 77^®) where 
the inverse abbreviates the solution of 

^ {E 6ss, - (f)) (F, F'; E) = S{f - F') <5,,, . 

Si 

(A2) 

The summation is over spin components, s = {t,)-}- For 
a non-local Hamiltonian an integration over the interme¬ 
diate real-space coordinates would also be present, and 
the remaining derivation is unaffected. 

The z-component of the spin density is given by 


TOz(F) = dE (G^^{r,f]E) - Gii{f,r;E)^ 


(A3) 

so we must solve for the diagonal (in spatial coordinates 
and spin labels) parts of the KS GF. 

Expanding Eq. (A2) we find the following 
among the spin blocks: 

relations 

{E - nffir)) Gn(F, F'; E) = S{f - F') 


+ nf^{r)G^^{f,r';E) 

, (A4a) 

{E -Hf^{f))GMr,f';E) = 0 


+ 77fi®(F)Gn(F,F';£;) 

, (A4b) 

{E - Gu(F, F'; E) = S{f - F') 


+ nf^if)Gnir,f';E) 

, (A4c) 

{E -H^f{f))Guir,f';E) = 0 


+ nf^{r)Gu{ry-,E) 

. (A4d) 

Defining auxiliary GF blocks, G^(F) = (A - 
and Gi{E) = (A — solutions of 

I 

{E - Hf^{f)) Gt(F, F'; E) = 5(F - F') 

, (A5a) 

{E - G^(F, F'; E) = <5(F - f') 

, (A5b) 


the previous set of equations can be rewritten as follows: 
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G'n(^, r'; E) = G^(r , r"; E) 

+ J dfi G^(f, ri;E) (ri) r E) , (A6a) 

Gi^{f,f';E) = 0 

+ J dfi Gi(f, fi-E) G^^(fi, f'; E) , (A6b) 

G^^(f,r';A) = G^(f,f';A) 
+JdfiGi{f,ri-E)Hff{fi)G^l{ri,f';E) , (A6c) 

Gn{r,f';E) = 0 

+JdriG^if,ri;E)Hf^{fi)Gii{fi,f';E) . (A6d) 

After using Eq. (A6b) and Eq. (A6d) in Eq. (A6a) and 
Eq. (A6c), respectively, we obtain 

Gtt(^, r"; E) = G^(f, f E) 

+ ^^2 Gt(r , ri; A) G^(ri, fz; E) 

x'Hff{f 2 )G^^{r 2 ,r'-,E) , (A7a) 

Gu{r,r'-,E) = G:,{r,f';E) 

(^i) Gt{ri,r2; E) 

xH^f{r2)G^^{r2,r'-E) , (A7b) 

and the solution of these two equations can be repre¬ 
sented in an abbreviated form as 

Gn(A) = {e- cm > (A8a) 

Gu(A) = {e- Gt(A) ■' . (A8b) 

Here the inverse operation is for the (r, r') dependence, 
as in Eq. (A5); the spin dependence was already solved 
for. 

By the following operator identity 

^-1 _ ^-1 = A-^{B - A) = B-\B - A) A-^ , 

(A9) 

the difference between the two GE blocks of interest is 
given by 

Gtt(r,r';A)-G^^(f,r';A) 

= ydfi Jdr 2 G^^(f,rq; A) A(fi,r 2 ; A) G;^(r 2 ,r';E) 

= ydri jdr2Gii{r,fi;E)A{ri,r2]E)G^^{f2,r';E) , 

(AlO) 

with the energy dependent splitting defined by 
A(f, r E) = (nf^^r) - )) S{r - r') 

+ nf^{f)G^{r,f';E)nf^^{r') 

- nffir) G^{f, r'-E) • (AH) 


Now the concrete form of the splitting is given, making 
use of the terms in Eq. (Al). This is especially mean¬ 
ingful if the xc field points along the z direction and is 
responsible for most of the spin splitting: 

A(r , r E) = 2 (^B^''{r)+^{r) 6{r- f') 

+ (C(r)L_)) G^(r , rA) (c(r')L++ ')) 

- (e(r)L++i?-‘(f))G^(f,rA)(e(r')T-+B“‘(r-")) 

(A12) 

with the combinations L± = iLy and B'^^ = H™* ± 
iBy^^. Defining the spin flip KS susceptibilities as 

xf- = i {xS - < + ix'l + Xj^*) , (A13) 

xti = i (xS + < - ix“ + Xj^*) , (A14) 

from the definition in terms of GEs and Pauli matrices 


X^I{u;) = -lfdETr 


aaG{E + u! + iO) Im G{E) 


+ aalmG{E) ap G{E — uj — iO) , 

(A15) 


we arrive at the magnetization sum rule, 


mz{r) = 2 ydr"x+® (r ,d";0) B^''{r') + Amz{r) 

= 2jdr'X-l{rA';0)B^%r') + Am,{r) , 

(A16) 

where the two equivalent variants follow from the two 
forms of the operator identity, Eq. (A9), and the small 
correction Amz(r) arises from the spin-orbit coupling 
and external magnetic field contributions to A(r , r E), 
{i.e. excluding the xc part). 


Appendix B: Kohn-Sham susceptibility in a basis 
derived from KKR scattering solutions 

In the atomic sphere approximation (ASA), the KKR 
GF has the form^^: 

Gt^if,r';E)=J2YL{f)Gti^^^^,ir,r';E)YL-ir') (Bl) 

LL' 


with 

GI^^^Ae r'-, E) = RUr-E) Gl^^^^E) R^,ir'; E) 

+ 5,,5LL'RUr^-E)Ht,{r^-E) . (B2) 

The position arguments of the GF are measured from 
the nearest atomic site, labelled i. The orientation is de¬ 
noted f and the length r; furthermore, r< = min{r, r'} 
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and r> = max{r, r'}. The angular dependence is ex¬ 
panded in spherical harmonics YL{r), with composite in¬ 
dex L = (£, m), and the two spin components are labelled 
by s. The radial functions E) and E) are so¬ 

lutions of the Schrodinger equation for the KS potential 
centered in site i for a given energy E, which are reg¬ 
ular and irregular at the nuclear position, respectively. 
The structural GF, jj^,{E), contains the information 
about the geometrical arrangement of the atomic sites 
and the multiple scattering contributions. We assumed 
a collinear magnetic state and no SOC, for simplicity of 
presentation. 

We use a basis of radial functions constructed from nor¬ 
malized regular scattering solutions computed at several 
energies Ei, within the range of the valence states (four 
energy values are sufficient for the basis construction): 




RUEEb) 

/o^drr2 Rl^{r-,EbY 


(B3) 


For fixed i and by orthogonalizing the overlap matrix 

Gw=l drr'^ (B4) 

Jo 

and keeping only the two largest eigenvalues, we form two 
linear combinations of the reference basis functions by 
using the respective eigenvectors, which become the basis 
functions for atom i and angular momentum channel £. 

The most general form of the KS GF (with SOG and 
non-collinear magnetism), in the KKR representation 
and in our chosen basis, is thus 

E) ='^^4’ub{''’)Glib^jL'b'{E)4>]iib'{R) ■ 

bV 

(B5) 

The KS susceptibility, Eq. (4), is then also naturally 
expressed in this basis, 






T '; wj = 2^ Yl, {f)YL^ {fWl.bS.rWt^bS-'^) 117^2, oL^L^b^b, ^%byWjl^by)YL2 if')yu if') 

Lisibi--- 


(B6) 


The size of the matrices involved is as follows. For the 
GF, the number of rows or columns is Na x x Ng x 
Nb = 1 X + y X 2 X 2 = 6A when taking only an 

adatom into account (Na = 1), with £max = 3 and two 
radial basis functions. For the KS susceptibility matrix, 
the number of rows or columns is Na x Na x (Nl x NgX 

NbY = 4 X 64^ = 16384, as we have four Pauli matrices 
[a = x,y,z,0), with the same assumptions as for the 
previous example. 


Appendix C: Dynamical magnetic susceptibility 
from the Landau-Lifshitz-Gilbert equation 


The aim is to linearize the LLG equation. 


dm ^ -;„fr m dm 

= -jrnx + r] — x— Cl 

at m at 


with the effective field 

i?®® = - ^ -I- bx {t) Cx + by (t) iy (C2) 

where is the static part, and bx{t) + by{t) Cy is 

the small time-dependent transverse part. Under these 


assumptions, the same applies to the magnetization, 


rh{t) =mxit)ex+my{t)ey +M^ez , (C3) 


with rux, TUy M. 

After Fourier transforming, ^ —>■ —ioj, the dynamical 
transverse magnetic susceptibility can be extracted from 
the LLG equation as 

, X _ MzOJo (1 + -bo + irjuj _ yi?®® 

x+- (^) - • 

(C4) 

This form can be used to fit the hrst-principles TDDFT 
data and extract all the parameters. 

The resonance peak location is obtained from 


y-Imx+- = 0 
da; 


(w — wo)^ + (??wo)^ — 2 u}{lo — Wo) = 0 
yR®® 


-^ ^max — /-- - ? 

which is the result quoted in the main text. 


(C5) 
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